Selective sweeps in growing microbial colonies 



Kirill S Korolev^'fll, Melanie J I Miiller^ ^ i^l, Nilay Karahan^, 
Andrew W Murray Oskar Hallatschek^ and David R 

Nelson^'2'3 

^ FAS Center for Systems Biofogy, Harvard University, Cambridge, Massachusetts 
02138, USA 

^ Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA 
^ Department of Molecular and Cellular Biology, Harvard University, Cambridge, 
Massachusetts 02138, USA 

^ Max Planck Research Group for Biological Physics and Evolutionary Dynamics, 
Max Planck Institute for Dynamics & Self-Organization (MPIDS), Gottingen, 
Germany 

E-mail: papers.korolev@gmail.com (K S Korolev), 

mmueller@physics.harvard.edu (M J 1 Miiller) , nelson@physics .harvard. edu 
(D R Nelson) 



Abstract. Evolutionary experiments with microbes are a powerful tool to study 
mutations and natural selection. These experiments, however, are often limited to 
the well-mixed environments of a test tube or a chemostat. Since spatial organization 
can significantly affect evolutionary dynamics, the need is growing for evolutionary 
experiments in spatially structured environments. The surface of a Petri dish provides 
such an environment, but a more detailed understanding of microbial growth on Petri 
dishes is necessary to interpret such experiments. We formulate a simple deterministic 
reaction-diffusion model, which successfully predicts the spatial patterns created by two 
competing species during colony expansion. We also derive the shape of these patterns 
analytically without relying on microscopic details of the model. In particular, we find 
that the relative fitness of two microbial strains can be estimated from the logarithmic 
spirals created by selective sweeps. The theory is tested with strains of the budding 
yeast Saccharomyces cerevisiae, for spatial competitions with different initial conditions 
and for a range of relative fitnesses. The reaction-diffusion model also connects the 
microscopic parameters like growth rates and diffusion constants with macroscopic 
spatial patterns and predicts the relationship between fitness in liquid cultures and 
on Petri dishes, which we confirmed experimentally. Spatial sector patterns therefore 
provide an alternative fitness assay to the commonly used liquid culture fitness assays. 
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1. Introduction 

Traditionally, the theory of evolution has been developed by analyzing phenotypes and 
genotypes found in natural populations and fossil records [H [2]. However, due to 
recent developments in microbiology and modern genetics, evolutionary experiments 
are becoming a valuable research tool [3j. Laboratory experiments hold great promise 
for uncovering basic evolutionary mechanisms by allowing us to observe evolution 
over time. More important, experiments, unlike evolution in natural populations, 
can be repeated systematically to distinguish between general principles and historical 
accidents. Microbes are particularly suited for evolutionary studies because they are 
relatively simple, reproduce and evolve rapidly, and can be easily modified using genetic 
engineering. Experiments with microorganisms could also provide insights into tumor 
growth, the spread of antibiotic resistance, and directed evolution of microbes to produce 
medicines or biofuels [H [5] . 

One potential drawback of evolutionary experiments is that they are conducted in 
artificial laboratory environments, which are quite different from the natural ecology 
of the species studied. The choice of the laboratory environment is therefore very 
important because it could affect both the nature of observed adaptations and the 
evolutionary dynamics. Most microbial experiments, of interest to us here, are 
conducted in the well-mixed environments of a chemostat or a test tube. These well- 
controlled environments allow researchers to compare experimental results to theoretical 
predictions, but it is important to ensure that such results are generic, not environment 
specific. Spatial structure, absent in well-mixed cultures, can significantly affect 
evolutionary dynamics O [TJ 13 El [10] ; therefore, it is important to carry out experiments 
with growth conditions that allow spatial inhomogeneities to form. The surface of a Petri 
dish is an easy-to-use environment, chemically similar to the environments of a test 
tube or a chemostat, yet capable of sustaining and preserving spatial structure during 
colony growth. Recently, several studies have used microorganisms in Petri dishes to 
study spatial patterning, mutations, and evolution [HI [121 [13 [H [HI [iHl M [El [HI 
[T8l [T9l [201 [211 [22] . Our focus here is on compact growth of colonies containing different 
genetic variants. For a comprehensive review emphasizing the beautiful dendritic growth 
patterns that can arise at low nutrient concentrations see Ref. [23] . 

An important obstacle to a wider use of spatial environments in evolutionary 
experiments is the limited theoretical understanding of how basic evolutionary processes 
play out in a spatial context. In particular, one must have a way to measure fitness 
to study evolution, but there are few models that relate microscopic parameters of the 
organisms to macroscopic quantities that can be easily measured in the laboratory. In 
this paper, we formulate a coarse-grained model of spatial competition that fills this gap 
and carry out microbial experiments to test the model's predictions. This model is based 
on deterministic reaction-diffusion equations, which describe short-range migration of 
the organisms and their competition. Our numerical results are further supported by a 
geometric argument, which does not rely on the detailed assumptions about microbial 
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growth and migration. Stochastic effects due to number ffuctuations (genetic drift) 
can be included in the model [151 1211 ISSj, but they do not play a central role in our 
analysis and are neglected for simplicity. For a study of spatial competition between two 
neutral bacterial strains, which is dominated by genetic drift, see Ref. [17]. Although 
we focus on competitions, the model is sufficiently general to describe mutualistic and 
antagonistic interactions as well. 

Compared to well-mixed populations, spatial populations have a wide range of 
initial conditions because one has to specify not only the relative fractions of genotypes 
but also their spatial distribution at time zero. We examine several experimentally 
interesting initial conditions and calculate how the spatial distributions of genotypes 
changes with time. We find that the shape of the resulting spatio-genetic patterns is 
determined by the expansion velocities ratio Vi/v2, where vi and V2 are the expansion 
velocities of isolated colonies composed exclusively of strain 1 or strain 2. Relative 
fitnesses of the genotypes can then be estimated by comparing experimentally observed 
spatio-genetic patterns to the theory. In particular, the selective advantage in this 
context can be defined as s = Vi/v2 — 1. As we show below, this definition is 
closely related to the traditional definition in terms of exponential growth rates. For 
linear inoculations, our results in the long time limit agree with Ref. [23], where a 
phenomenological model of the patterns was first proposed. In contrast to Ref. |24j . 
we also study circular inoculations and carry out experiments to test quantitative 
predictions of our model. 

The theory can be directly compared to the experiments because the spatial 
distribution of genotypes on a Petri dish can be visualized and quantified with fluorescent 
markers [15]. Here, we briefly describe this technique; see the supplementary information 
(section SI) for more details. Microbial strains of interest are genetically modified to 
constitutively produce a fluorescent protein. The emission spectra of the proteins must 
be sufficiently different in order to distinguish the strains on a Petri dish. To study 
selective sweeps, a mixture of the strains is prepared in liquid medium, usually with 
the fitter strain in the minority. A drop of this mixture is then deposited on a small 
region of a Petri dish with solid growth medium. Different shapes of this drop lead to 
different initial conditions. Circular drops are naturally created by the surface tension 
forces when small drops of fluid are placed on a Petri dish. Linear drops can be created 
by gently touching the surface of the medium with a razor blade after dipping it in the 
mixture of the strains. These drops dry quickly, and microbial colonies start to grow, 
expanding by about a centimeter a week. The spatial distribution of genotypes can be 
observed during this expansion by fluorescent microscopy, as shown in figure [H 

Sectors in figure [1] are at the center of this study. In our experiments, cells are 
nonmotile and grow primarily at the expanding frontier; behind the front, the growth 
is limited by the lack of nutrients [271 EH]- As a result, the genetic composition in the 
interior of the colony does not change with time. The spatio-genetic pattern shown in 
figure m is then a frozen record of the temporal changes in the spatial distribution of 
genotypes at the expanding frontier. In figure [1^, we show growth of two strains with 
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Figure 1. (Colour online) Comparison of spatial segregation during a range expansion 
of Baker's yeast, Saccharomyces cerevisiae with (a) equal and (b) different growth 
rates of the two competing strains. The Petri dishes were inoeulated with a well- 
mixed population occupying a narrow horizontal linear region at the bottom of the 
images from which the sectors appear. As the populations expand, they segregate 
into well defined domains. Different colours label different genotypes. In (a), the two 
strains (yellow and blue) have the same fitness and the demixing is driven primarily 
by number fluctuations (genetic drift) [15l[6]. It is likely that the small variations with 
horizontal position in boundary slopes are related to undulations of linear fronts, which 
are hard to suppress when the front is very long [55]. In (b), the sector is formed by 
the fitter strain (black) , and the sector expansion is caused by the difference in growth 
rates of the strains, or, in other words, by natural selection. In both (a) and (b), the 
scale bars are 500 /im. 



(a) 


(b) 


(c) 


(d) 







Figure 2. (Colour online) Colony edge at single cell resolution (mature yeast cells are 
5 /Ltm in diameter), (a), (b), (c), and (d) are successive images (at two hour intervals) 
of the same region near the edge of a growing S. cerevisiae colony inoculated with a 
razor blade. Note the formation of a green (light gray) sector on the lower right. The 
two strains have approximately the same fitness in this experiment. 

the same fitness that differ only in the colour of a fluorescent marker. Before deposition 
on the Petri dish, these strains were combined in a 1:50 ratio and thoroughly mixed. 
However, the two different colours (genotypes) demix and form sectors. This demixing 
is caused by number fluctuations (genetic drift) at the expanding edge [IHl E]. One 
can see this stochastic process at the resolution of a single cell (5 /im in diameter) 
in figure |2l In contrast to neutral demixing, figure [T|3 shows sector formation in a 
colony founded by two strains with different fitnesses; the fitter strain is in the minority 
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initially. Over time, the fitter strain displaces the other strain, increasing its share of 
the expanding frontier. This expansion is also subject to number fluctuations, which 
are responsible for the sector boundary wiggles, but the average shape of the sector is 
determined primarily by the deterministic force of natural selection. For the purpose 
of measuring relative fitness, the effects of genetic drift can be averaged out, given a 
sufficient number of repeated experiments. Our deterministic model strives to describe 
this average sector shape and is not capable of describing the randomness in boundary 
motion or the dynamics when the number of the fitter organisms is so small that the 
number fluctuations can lead to their extinction, an effect discussed in Ref. [23]. 

This paper is organized as follows. We formulate a competition model assuming a 
well-mixed environment, such as a mechanically shaken test tube, in Sec. [21 Under 
these conditions, a particular microbe visits virtually every region of the carrier 
fluid in a cell division time, and the system is effectively "zero-dimensional." This 
model is then extended to account for lateral migrations during a range expansion 
in Sec. [3] In Sec. [31 we analytically derive the spatial patterns created by two- 
species competition (including results for colliding circular colonies) using a very general 
argument, which does not rely on the microscopic details of microbial growth and 
migration. The theoretical predictions are then compared to experiments in Sec. [5] 
Concluding remarks are contained in Sec. [6] The supplementary information (section SI) 
contains the experimental and numerical methods as well as additional data supporting 
our conclusions. 

2. Modeling competition in a well-mixed environment 

A competition experiment is a standard way to measure relative fitness of two microbial 
strains in a well- mixed environment. During a competition experiment, the stains are 
introduced into a fresh medium, and their relative abundance is measured over time. 
Initially, the number of cells grows exponentially, but the growth eventually slows down 
as the system approaches the stationary phase due to crowded conditions. This behavior 
is captured by a simple Lotka-Volterra-type model [29] : 




which is the most general model with quadratic nonlinearities. Here, ci{t) and C2(t) 
are the concentrations (number of cells per unit volume) of strain 1 and strain 2 
respectively. The constants gi and g2 are their exponential growth rates; and the 
constant matrix dij > describes nonlinear interactions. For a mono-culture consisting 
of a single yeast strain, Refs. [301 [SI] showed that population growth can be described 
by the logistic equation (i.e. with quadratic nonlinearities) very accurately. 
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Figure 3. Schematic representation of four different behaviors of equation ^ in phase 
space, (a) ei > and €2 > 0. (b) ei > and £2 < 0. (c) ei < and ti > 0. (d) ei < 
and 62 < 0. The arrows represent the direction of trajectories in the phase space, and 
the dots represent the fixed points. 



By rescaling c^'^^it) = diiCi{t)/gi and c^'^^it) = d22C2{t)/g2, we can recast 
equation ([1]) in a slightly more convenient form 

^ci(t) = gMt)[l - ci{t) - C2{t)] + eici(t)c2(t), 
d 

-C2{t) = g2C2{t)[l - Ci{t) - C2{t)] + e2C2(t)Ci(t), 

dt 

where ei = gi — (ii2fi'2/'^22 and €2 = g2 — c?2i5'i/'^ii5 and we use the same symbols Ci{t) 
and C2(t) for the rescaled concentrations. When ci(t) + C2(t) ^ 1, we shall say that 
the competition takes place "under crowded conditions" . Our notation emphasizes that 
the initial stage of exponential growth is usually much shorter than the second phase 
of competition under crowded conditions, i.e. the growth rates gi > and g2 > 0, are 
typically much larger than the small quantities |ei| and |e2|. When > 0, strain i grows 
faster in the presence of the other strain, e.g. by feeding off an excess production of 
a useful amino acid. When < 0, strain i grows slower in the presence of the other 
strain, e.g., because of a secreted poison. 

Note that equations ([2]) always have at least three fixed points: (0,0), (1,0), 
and (0,1). In addition, for some values of the parameters, there is another fixed 
point (cJ,C2) in the physically relevant domain of nonnegative Ci and €2- The fixed 
point at the origin is always unstable. 

Five different behaviors are possible depending on the values of ei and €2] four of 
these are illustrated in figure [31 If ei and €2 are positive, the interaction of the species 
is mutualistic (i.e. the presence of the first strain helps the second strain grow and vice 
versa) and leads to a single stable fixed point (0^,03) with a nonzero concentration of 
both strains. Since the focus of this paper is on competition, we do not pursue this 
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Figure 4. (Colour online) Solutions of equation ([2]) plotted in phase space (ci,C2) for 
negative ei and positive 62- The arrows indicate the direction of time. For the initial 
conditions, we chose positive ci(0) and 02(0) close to zero with each trajectory having 
a different value of ci(0)/c2(0). In this plot, gi — 1.5, 52 = li = —0.1, 62 = 0.6. 
Note that, initially, all trajectories bend toward the ci-axis (increasing C1/C2), but the 
first strain is nevertheless eliminated at long times, as the system approaches the stable 
fixed point (0, 1). 

cooperative possibility further here. If ei and €2 are negative, there are two stable fixed 
points: one with Ci = 1 and C2 = 0, the other with Ci = and C2 = 1. The system 
reaches one of these fixed points depending on which strain is more prevalent initially. 
An incoming separatrix divides the phase space in two domains of attraction and feeds 
into an unstable fixed point with > and C2 > 0. If ei > and €2 < 0, there is only 
one stable fixed point ci = 1 and C2 = 0. Similarly, ci = and C2 = 1 is the only stable 
fixed point when ei < and 62 > 0. Finally, there is a degenerate case ei = €2 = when 
the dynamics is determined only by the logistical growth. In this case, depending on 
the initial conditions, the system lands somewhere along the line of neutral fixed points 
defined by Ci + C2 = 1 . 

Interestingly, when 6162 < 0, the ultimate result of the competition is independent 
of the exponential growth rates gi and g2. The dynamical path, however, does depend 
on these parameters. For example, for ei < < 62 and gi > g2 (the situation of 
particular interest to us here), Ci(t) initially increases, but eventually falls off as the 
total population density rises settling at (01,02) = (0,1), see figure HI This behavior 
has important implications for the competition at the frontier of a colony expanding 
in space, where the concentration of cells is always low. We show in Sec. |3] that the 
outcome of a spatial competition experiment is determined primarily by the exponential 
growth rates gi and g2 rather than ei and 62, in contrast to the aforementioned behavior 
in a well-mixed population. 
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3. Competition and migration 

In this section, the competition model formulated in Sec. [2] is generalized to spatially 
extended environments. To this end, we need a model of cell migration. In general, cell 
migration is a complicated process, which could involve chemotaxi^, swarming, and 
random wandering. Although all of these can be important biologically, some can be 
neglected in appropriately designed spatial competition experiments. For example, cells 
of Baker's yeast, Saccharomyces cerevisiae, often used in microbiological experiments, 
are nonmotile, and many bacterial cells, e.g. Escherichia coli, can be prevented from 
swimming by eliminating functioning flagella or using a high concentration of the 
gelation agent in the growth medium. 

For nonmotile microbes, the only mechanism of cell migration is cells' pushing on 
each other as they increase in size before cell division. Even a cell division of an isolated 
cell leads to migration because at least one of the offspring is generally displaced in 
a random direction relative to the position of its parent. At the colony front, where 
cells are relatively free to move, cell migration can be approximated by a random-walk- 
type process caused by growing cells pushing each other in random directions. The 
diffusion constant of such random walks must depend on the local concentration of cells 
because sector boundaries in the interior of a colony like the ones shown in figure [1] 
do not change with time [15]. For a concentration-independent diffusion constant, 
the boundaries would slowly disappear as cells of different colours gradually mix [2^ : 
therefore, migration must be arrested at high cell concentrations. A possible mechanism 
of this arrest is a significant reduction in growth rate due to nutrient depletion. At the 
edge of the colony, however, the cell density is low, and cells move readily due to the 
jostling caused by cell growth, as is evident from the wandering of the sector boundaries. 

Since the exact dependence of the spatial diffusion constant on the local cell 
concentration is unknown, we have explored a family of functional forms, namely 



where -Dqi is the diffusion constant of the first strain in the limit of small local 
cell concentration, and cti is an adjustable parameter that allows us to explore the 
sensitivity to rapid variations in the rate of migration near the frontier. Here, Ci{t,x) 
and C2{t,x) are cell densities per unit area, rescaled as in Sec. El An analogous 
dependence, -D2(ci, C2) = -Do2(l — ci — C2)°^6'(l — Ci — C2), is assumed for the second strain; 
here 6{x) is the Heaviside step function, 6{x) = 1, x > 0, and 6{x) = otherwise. The 
choice of a monotonic dependence of the diffusion constants on cell density is motivated 
by the monotonically decreasing supply of nutrients from the outside to the inside of 
the colony. To check whether a non-monotonic dependence would affect our results, we 
also explored -Di(ci, C2) = -D2(ci, C2) = -Do(ci + C2)(l — Ci — C2)6'(l — ci — C2). Although 

II Chemotaxis is an ability of cells to direct their motion in response to a chemical signal, e.g. food or 



^i(ci, C2) = Doi{l - Ci - C2)"' for ci + C2 < 1, 
-Di(ci, C2) = for ci + C2 > 1, 



(3) 



toxins. 
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Figure 5. (Colour online) Both (a) and (b) show the solution of equation Q plotted 
for gi = 11.8, 32 = 10, ei = ea = 0, ai = a2 = 1, Dqi = 1.18 • lO'^ and £'02 = 10"^ 
note that we chose gi/g2 = -D01/D02 = 1-18 to facilitate comparison with figurc[T}D (see 
discussion in text). In these units, the habitat is a 1 x 2 rectangle and is initially empty; 
only the bottom sixty percent of the habitat is shown because the top part remains 
empty throughout the expansion. The origin of the expansion is a line at the bottom 
edge of the images, where we impose the boundary conditions that C2{t,x,0) = 1 
and ci(t,x,0) = except in the width 2^* region near the center of this boundary, 
where C2{t, x, 0) = and ci(t, a;, 0) = 1. No- flux boundary conditions are imposed along 
all other edges. Equation (|4]) is solved on a square grid of 256 x 512 points, (a) The 
concentration of the first strain is shown in red (dark gray) and of the second strain in 
green (light gray). The maximal colour intensity corresponds to the concentration of 1, 
and the lowest to the concentration of 0. This colour scheme is chosen to facilitate the 
comparison with the experimental data shown in figure [TJ (b) The same solution as in 
(a) , but only the concentration of the first strain is shown to highlight its establishment 
as a sector early in the expansion. Brighter regions correspond to higher concentration 
of the first (red) strain. 



the speed of population waves is different in this model, the shape of spatial patterns 
and their connection to relative fitness remain the same; see supplementary information 
(section S4). 

Our model of spatial competition then takes the following form 
' dci{t, x) 



dt 



dc2(t,x 



V ■[Di{ci,C2)Vci{t,x)] + 
gici{l - ci - C2) + eiCiC2, 

- = V - [D2iCi,C2)VC2{t,x)] + 
9202(1 - Ci - C2) + e2C2Ci. 



(4) 



The initial and boundary conditions are chosen to mimic selective sweeps in microbial 
colonies. For example, in figure 13 only a very small region is initially occupied by the 
advantageous genotype, which corresponds to the time when a sector begins to appear 
in figure [T|d. Our model cannot describe the earlier part of the range expansion when 
number fluctuations are important, because equations (jlj) are deterministic and treat 
the cell densities as continuous functions. During this early stage, the advantageous 
genotype becomes extinct stochastically everywhere but a few spatial locations, where 
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Figure 6. (Colour online) The solution of equation Q plotted for gi = 12.5, 52 = 
10, ei = £2 = 0, ai = q;2 = 1, and Dqi = ^02 = 0.01. The habitat is a 10 x 10 square 
inoculated with a circular drop of radius 2 shown in gray. We assign ci and C2 in the 
initial circular boundary to be either or 1 (in blocks) to mimic the sectoring pattern 
produced by a short period of genetic drift with a relatively small selective advantage. 
No-flux boundary conditions are imposed along all edges. Equation ([4]) is solved on a 
square grid of 2560 x 2560 points. The concentration of the first strain is shown in red 
and of the second strain in green. The maximal intensity (of red or green) corresponds 
to the concentration of 1, and the lowest to the concentration of 0. 

it gives rise to small sectors. These sectors can then be used as the initial condition 
in our deterministic model. We solve equations numerically; see the supplementary 
information (section SI) and figures [5] and [6l Note that our model of range expansions 
has two spatial dimensions, while, in experiments, colonies also gradually thicken in the 
direction perpendicular to the plate |32]. We neglect this gradual thickening here. 

From figures [lb and [5], one can see that there are at least two stages in the sector 
formation. During the late stage, the two sector boundaries are far apart, and the 
interior of the sector is dominated by the advantageous strain. At this stage, the sector 
boundaries maintain a constant angle with the direction of the expansion, which we 
explain in Sec. H] using a very general geometric argument. During the early stage, the 
size of the sector is comparable to the width of sector boundaries, and the two boundaries 
interact. By neglecting the nonlinear terms in equations (j4]) near the frontier, where Ci 
and C2 are small, we can qualitatively understand how the duration of the first stage 
depends on the relative fitness when the fitness difference is small (see figure [7]). In 
this case, we can neglect the bulging of the sector and assume that the population 
wave front is approximately flat. For simplicity, we also assume Dqi = Dq2 = Dq. 
The X-axis is taken to be along the front and the y-axis to be perpendicular to the 
front. For the Fisher equation [33], a steady state is reached in a frame co-moving 
with the population wave {x' = x and y' = y — 2y/g2Dot) when only the second strain 
is present. In this reference frame, (72C2 ~ —Dq^^ — (note that the nonlinear 

terms are neglected). Since the fltness difference and the concentration of the flrst 
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Figure 7. (Colour online) The length of the initial stage of interacting sector 
boundaries from equation as a function of the ratio of expansion velocities Vi/v2- 
The quantity yi^2 is the distance from the origin of the sector to the closest point 
where ci = 1/2. Here, we vary gi while keeping Dq = 10^^ and 52 = 10 fixed. For 
small fitness differences, we expect 2/1/2 ~ (51 — 52)^^ from equation The data 
from numerical solutions of equation Q is shown as dots, and the solid line is a fit 
to A/ (gi — (72) + B, where A and B are fitting parameters. 



strain are small, the dynamics of the second strain along the y-axis is approximately 
unchanged. Moreover, since the spatial distribution of the first strain along the y- 
axis is the same as that of the second strain, a similar equality holds for the first 
strain g2Ci ^ — Dof^ ~ ^2^- Upon using these two observations and the linearized 
version of equations (jl]), we find that the dynamics of the first strain near the frontier 
is approximately given by a linear diffusion equation with a source, 

Ol = ^0 9^2 + i9i-92)ci{t,x,y), (5) 

Note that this equation is invariant with respect to translations along y', so we can 
treat the frontier as a quasi-one- dimensional population. For ci(0,x',y') = S{x'), 
corresponding to a point-like inoculant at the frontier, the solution of equation ([5]) 
is 

rJi^^') = — ^ p(gi-g2)t^-4i^ /g) 

Therefore, the characteristic time necessary for the first strain to dominate the sector 
scales as {gi — g2)~^- This divergent time scale is indeed observed in the numerical 
solutions; see figure [71 

The dependence of spatial patterns during two-species competition on various 
parameters in equation (jlj) was investigated in the context of linear expansions. We 
varied the exponents ai and ^2 by factors of 2, ai = 0^2 G {1/4,1/2,1,2,4}; the 
diffusion constants Dqi and Dq2 by factors of 10, Dqi = Dq2 G {10~^, 10"'^, 10~^}; 
and the growth rates gi and g2 by factors of 10, gi G {1,10,10^} with gi/g2 G 
{1, 1.01, 1.05, 1.1, 1.3, 1.5, 2}. The competition parameters ei and €2 were independently 
varied relative to minjgfi, (72}, ei,e2 G [— 0.5, 0.5]min{(7i, (72}. This numeric exploration 
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helped us identify important parameter combinations that control the shape of spatial 
patterns. We now turn to the discussion of these results. 

From numerical analysis, we made an important observation that the expansion 
velocity of a strain growing in the absence of the other strain depends only on the 
exponential growth rate and diffusion constant right at the frontier and is given by 

vi = 2^/Doigi, 

independent of ei, €2, ai, a2 (see supplementary information section S4) and in 
agreement with the classic Fisher-Kolmogorov wave theory [331 IM], which provides 
an exact solution for a simpler model with concentration-independent diffusivity. This 
agreement is not surprising because the speed of a Fisher population wave is determined 
only by the dynamics at the foot of the wave front [331129]. and equation (jHl) ensures that 
the diffusivity approaches a constant for small ci, C2. The intuition behind equation d?]) 
is that the wave speed depends both on the growth rate [gi) and on the rate of undirected 
migration (-Dqi) that brings cells to unoccupied territories. The detailed shape of 
the wave front, however, does depend on parameters such as ai and gi/D^i (or 02 
and 5'2/-Do2)- The parameters ei and €2 are irrelevant when only one strain is present 
because C1C2 = in this case. 

Even when both strains are present, the knowledge of expansion velocities vi and V2 
is sufficient to describe the major (large scale) features of the resulting spatio-genetic 
pattern. To see this, note that the behavior of population fronts far from the sector 
boundaries is the same as when only one strain is present because the concentration 
of the other strain vanishes away from the boundaries. We also found that the initial 
position of the sector boundaries is determined only by vi and V2 (and independent of ei 
and €2) because, at the tip of the advancing front, the product C1C2 is exponentially small 
compared to Ci and 02- On smaller length scales of the order ■^/Dqi/S'i (or \/ 92)-, all 
parameters play a role. In particular, all parameters affect the shape of concentration 
profiles and the position of sector boundaries during the early stage of sector formation. 

Note that, with nonzero ei and €2, sector boundaries still move behind the wave 
front even though the diffusivity is zero. To understand this motion, consider strains 
with g\ > g2 and ei < < €2. The first strain expands faster, but, after the front has 
passed, the dynamics under crowded conditions, discussed in Sec. El favors the second 
strain. As a result, any region with a nonzero concentration of the second strain is 
eventually colonized by it. This behind-the- front competition should lead to a finite 
displacement of the boundary because every sector boundary has a finite width due to 
the discreteness of the number of organisms; see figure [HI Since, for many microbial 
strains, the boundary width is small, and (71, (72 ^ l^i I) 1^215 we do not expect to observe 
this type of sector boundary displacement (very different from a Fisher genetic wave) 
experimentally. 

What is the relation between Vi/v2 and relative fitness in liquid cultures? In well- 
mixed populations, selective advantage is often defined from the ratio of exponential 
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Figure 8. Schematic illustration of transverse boundary motion behind a front 
advancing in the y-direction for gi > g2 and ei < < £2- The plots show concentration 
profiles along a linear cut along the x-direction, parallel to the front and inoculant. (a) 
The concentration profiles a short distance behind the population frontier. There is an 
overlap region, where both ci{x) (dashed line) and C2{x) (solid line) are nonzero. This 
region has a finite width (from A to B) because the discreetness of the number of cells 
is inconsistent with infinitesimally small values of the concentrations. (Sec Ref. |35j 
for a more detailed discussion of this issue in a related model.) (b) The concentration 
profile at the same spatial location as in (a), but after a very long time. The interval 
between A and B is now occupied exclusively by the cells of the second strain, which 
wins out under crowded conditions. 



growth rates, = 9i/ 92 — ^- At expanding frontiers, we define s = Vi/v2 — 1 by 
analogy. There may not be a direct correspondence between s and Swm because the 
former involves diffusion constants Dqi and -D02 in addition to the growth rates gi 
and g2- For example, a mutation providing a means of motility could be beneficial in a 
Petri dish (due to faster spreading) and deleterious in liquid (due to its metabolic cost). 
However, two special cases Dqi = Dq2 and Dqi ~ gi (with Dq2 ~ (72) are of interest. 

Equal diffusion constants should be a good approximation for mutations that affect 
growth rate, but do not affect motility directly, which is possible when motility does not 
depend on the growth rate strongly. Growth independent motility was, e.g., observed 
in swimming Bacillus subtilis cells [H]. Under these assumptions, 1 + s„m = (1 + s)^, 
and = 2s for s <^ 1; see equation ([7]). The other possibility Doi ~ g^ could be 
a good approximation when motility and growth are strongly linked. For example, 
colonies of S. cerevisiae studied here expand due to cell growth, and it is reasonable to 
assume that DQi ~ g^a , where a is the average cell size. In this follows 
from equation (I7j) ; see figure [9l 

Our simple definition of s (s = fi/f2 — 1) has three advantages: expansion velocities 
can be easily measured, larger expansion velocity results in greater colonized territory 
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Figure 9. (Colour online) The comparison between selective advantage in liquid 
culture and on Petri dishes within the reaction-diffusion model. The red circles show 
the results of the numerical solution of equation ^ for single strain expansions. The 
black line shows the theoretically predicted linear dependence. We varied gi to mimic 
different growth rates in the experiments and used 52 = 10,ei = e2=0,Q:i=Q!2 = 
1, Dqi = 10~^gi, and Z)o2 = 10~^. In these units, the habitat was a 1 x 10 rectangle 
and was initially empty. Each expansion was started at the shorter edge of the habitat, 
where we imposed the Dirichlet boundary condition forcing strain density to be 1. No- 
flux (Neumann) boundary conditions were imposed along all other edges. The grid 
size used was 128 x 1280 points. 

and, therefore, greater access to nutrients, and the ratio Vi/v2 is closely related to the 
traditional definition of obtained from the exponential phase of well-mixed cultures. 

In summary, the main conclusion of the mechanistic modeling embodied by 
equation (jlj) is that the essential features of spatio-genetic patterns formed during range 
expansion are insensitive to the details of the model and are determined by a single 
dimensionless parameter Vi/v2 (the ratio of expansion velocities). We further support 
this conclusion in the next section by deriving the shapes of the spatial pattern without 
relying on the microscopic dynamics of growth and migration. Additional tests of the 
robustness of our reaction-diffusion model to changes in the modeling assumptions (such 
as varying the form of the concentration-dependent diffusion constants) are presented 
in section S4 of the supplementary materials. 

4. Sector shapes and the equal-time argument 

In this section, we develop an analytic argument to understand selective sweeps and 
explore different ways of measuring relative fitness from macroscopic competition 
experiments. This argument relies on two assumptions: the strains expand with constant 
velocities and a patch occupied by one strain is impenetrable to the other strain. We 
primarily focus on two experimental geometries: linear and circular. 
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Figure 10. (Colour online) Equal-time argument and sector shape in a linear 
geometry. The wiggles in the sector boundaries represent genetic drift, neglected in 
most of this paper. 



4.1. Linear inoculations 

Upon establishment, i.e. when sector boundaries are sufficiently far apart, the sectors 
from a linear inoculation have a triangular shape. Following Ref. [21], we explain this 
shape by the simple equal-time argument illustrated in figure [TOl It should take the 
same amount of time for the second (less fit) strain to grow along the expansion direction 
as it takes the fitter first strain to grow a longer distance along the boundary. More 
generally, for any point along the front, we can define a length p that is the length of the 
shortest path connecting this point and the origin of the expansion and lying entirely 
in the territory occupied by the same strain. (This path is a straight line for linear 
inoculations, but curved for radial ones, as we show below.) Then, the ratio of p and 
the appropriate expansion velocity is the time necessary to form this particular spatial 
pattern. This time is the age of the colony since inoculation and must be the same for 
all points along the front. From this observation (see figure [TOj) . we conclude that the 
bulging shape of an advantageous sector is an arc of a circle of radius Vit and angle 
given by 



which is equivalent to the result obtained in Ref. |23]. 

Our arguments are somewhat analogous to the Huygens-Fresnel principle in 
optics [36] • Each point along the front has a potential to create an outgoing circular 
population wave that spreads with velocity Vi or V2, and unoccupied territories are 
colonized by the strain that gets there ffist. Unlike in optics, this principle can only 
be used to construct wave fronts at infinitesimal time steps because, once a region is 
colonized by one strain, it becomes impenetrable to the other strain. 

The equal-time argument breaks down on length scales smaller than Dq/vi oc 
^JDq/ gi because the Fisher velocity will in general depend on the curvature of the 
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Figure 11. (Colour online) Equal-time argument and sector shape for a circular 
inoculant of radius Rq. 



wave front. For example, advancing a circular front requires more time than advancing 
a linear front because more cell divisions are necessary to cover the larger area colonized 
by the curved front. The Fisher velocity approaches the limiting value, equation ([7]), 
provided the local radius of curvature is much larger than Dq/vi p9]. Dq/vi is also the 
characteristic width of a sector boundary, where the equal-time argument breaks down 
because the strains are intermixed. 



4-. 2. Circular expansions 

We now turn to radial range expansions, resulting from moderate-sized circular pioneer 
populations, e.g., created by placing a drop of the inoculant on the surface of a Petri 
dish. Due to surface tension, small drops (of order 5 /zl to 3 mm in diameter) are 
naturally circular. A precisely defined initial shape is a significant advantage over linear 
geometries susceptible to front undulations. This advantage is further strengthened 
because a two-dimensional reaction-diffusion population wave started from an irregular 
island of cells becomes more and more circular as the expansion continues. Even more 
important, the interactions of yeast cells lead to an effective surface tension suppressing 
front undulations |32] . 

The equal-time argument also yields the shape of sectors in the circular geometry; 
see the schematic plot (figure [11]) and numerical solution of equation (figure |6]). 
Equating the infinitesimal time increments along the radius of the wild-type colony and 
the curved sector boundary leads to a differential equation: 



dr y^dr2~+"(rd^p 

formulated in polar coordinates (r, 0) with the origin at the center of the expansion. 
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Figure 12. (Colour online) Illustration of the equal-time argument for the shape of 
the bulge induced by a faster growing strain (strain 1) in the circular geometry; see 
also figure [TTJ The central part of the bulge (between C and D) is an arc of a circle 
bounded by two tangents CK and DK (black dashed lines) to the sector boundaries 
at their origin. The rest of the bulge (between F and C, and between D and E) is 
described parametrically by equation (jl2p , where, for any point A on this part of the 
bulge, the parameter p is the distance between the center of the homeland and the 
intersection point B between the closest sector boundary and its tangent AB passing 
through A. Equating the total expansion time of the fitter strain, first, along the sector 
boundary KB and, then, along the tangent AB to the expansion time of the other 
strain along the radius of the green segment immediately yields equation ([12]) • 



where the different signs corresponds to boundaries turning cloclcwise and countercloclc- 
wise, and Rq is the initial radius of the population. Thus, we can identify sector bound- 
aries in microbiology with the famous logarithmic spiral of Bernoulli, which also de- 
scribes the Nautilus shell and insect flight patterns [37] . 

The shape of the bulge at the frontier can also be calculated with the equal-time 
argument. The top of the bulge, well away from the sector boundaries, is a circular arc 
of radius vit centered on the origin of the sector; the angular length (pt, of this region (see 
figure [12]) is given by 



just as in the linear geometry. Beyond </'b/2, the bulge is closer to the sector origin 
than vit because these points along the bulge cannot be connected to the origin of the 
sector by a straight line without intersecting the territories occupied by the other strain. 
In this case, the shortest allowed path back to the founding population is a straight line 
passing through a specific point (point B in figure [I2]) followed by a curved path along 



The solution reads 




(10) 
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the sector boundary to the origin of the sector. From figure [T^l one can easily see 
that the straight path is tangent to the sector boundary. Upon invoking the equal-time 
argument (see figure fT2l) . we find that the shape of the bulge beyond 0;,/2 is described 
by 



X = p sin K + (_Ro + — p) (sin k + 



y = P cos K + {Rq + V2t — p) (cos K — 



where t = is the time of the inoculation, n is given by 




1 cos k) 



1 sin k) 



(12) 



K 



'4-llnfA 

v'2 \Ro 



(13) 



and p G (-Rq, Rq + V2t) is a parameter equal to the length of OB in figure | 

Since the two sector boundaries turn in opposite directions, they must eventually 
meet, thus, enclosing the less fit strain within the population of a faster growing strain, 
as shown in figure [T^ This enclosure occurs at = ivr, and, from equation ([TU]) . we 
calculate the distance Rj from the center of the inoculation to the point where the two 
boundaries finally meet. 



Rf = Rq exp 



'n-V2 



Ro exp 




(14) 



where the last equality follows from our definition of selective advantage Wi = (1 + s)v2- 
Note that the time tf = Rf /v2 at which the second strain is enclosed by the first strain 
is exponentially large as s — )■ 0. Therefore, for small s, a competitive exclusion requires 
a much longer time than {gi — (72)^^ = 1/swm predicted by the well- mixed population 
models. 

Equation (fT2|) is valid only up to t = Rf/v2, and, after the enclosure, the shape of 
the colony is determined by the expansion of the first strain with velocity vi. 

Another interesting consequence of the equal-time argument is that the point where 
the two strains meet at the colony edge moves tangentially to the front of the less fit 
strain with a constant velocity v± = \/ v\ — v\\ see figures [10 and [HI Indeed, the 
distance /(t) between the two sector boundaries at a linear front increases with time 
as (equation [8]) 



2t'2(t-to)sin(0/2) 
2. 



l{t) = l{to 

= l{to) + 2^vl-vl{t-t^), 



(15) 



and the angular size of a sector if{t) in the circular geometry grows as (equation [TO 

2^- 



(p{t) = ip{to) 



In 



V2 



m 

R{to) 



(16) 
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(a) 



(b) 




(d) 



Figure 13. The enclosure of a slower growing strain (initially in the majority) by 
the faster growing strain during a competition experiment. The fitter strain is shown 
in gray, and the other strain is shown in black. Four consecutive snapshots of the 
numerical solution of equation (U) are shown in (a), (b), (c), and (d). The enclosure 
leading to the heart-shape occurs shortly before the snapshot shown in (d). 



where the factors of 2 are due to boundary motion at both edges of a sector. The lateral 
expansion velocity is then given by 



Thus, v± is a constant. One can also show that v± = Vi — V2 and these three vectors 
make a right triangle {v± _L V2); here Vi and V2 are front velocities at the sector 
boundary. Sectors of deleterious strain also obey equation f ll7p . but with a negative 



4-3. Colony collisions 

We conclude this section by considering competition between two strains not initially 
in contact. At the beginning of the experiment, two circular colonies inoculated with 
different strains grow independently each with their own velocity. Eventually, however, 
the colonies collide; see figure [TH For simplicity, we assume that the initial radius of 
the colonies is much smaller than the distance / between the colonies. Provided the 
nutrients remain abundant by time t, the colony boundaries are circles with radii Vit 
and V2t except for the collision boundary. The shape of the collision boundary follows 
from the equal-time argument; see figure UM Remarkably, this boundary is also a circle, 
with radius Rb given by 




(17) 



lateral velocity, v± = ~^/\vf 




Rb — I 
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Figure 14. (Colour online) Schematic picture of a colony collision. 



Note that Ri, diverges as the selective advantage s 0. The center of this circle is 
located on the line connecting the centers of the colonies, distance xo away from the 
center of the colony established by the first (faster growing) strain. We find 



where positive xq corresponds to the direction towards the colony with the second strain, 
and negative xq corresponds to the opposite direction (see figure UM . Similar to the 
selective sweep in the circular geometry, the less fit strain is eventually enclosed by the 
other strain. The time to this enclosure is given by 

= ^t±f^ = ±, (20) 

V2 V2S 

as one can see from figure [Ml and equations (ITSl) and (fT9l) . 

We derived sector and colony shapes using the equal-time argument and the results 
are consistent with our microscopic model. Because this argument relies only on the 
assumptions of constant expansion velocities and impenetrability of occupied regions, 
the equal-time argument applies more generally. As a result, we expect that the 
competition outcome is determined only by the ratio of expansion velocities of the 
two strains or species for any model consistent with these two assumptions. 

5. Comparison with experiments 

In the two preceding sections, we described theoretical predictions for the patterns 
of genetic diversity arising from spatial competitions. In section [3l we developed a 
generic reaction-diffusion model in order to predict macroscopic spatial patterns from 
microscopic parameters like cellular growth rates and effective diffusion constants. 
In section IH we then derived analytical formulae for boundaries between regions 
occupied by the competing strains using the equal-time argument. In this section, 
we experimentally test our theoretical predictions, using the expansion of the budding 
yeast Saccharomyces cerevisiae on agar surfaces as a model system. We also employ the 
analytical formulae from the equal-time argument, equations ([8]), (fTOl) . (fT8|) . and ( IT9|) . 
as a means to measure relative fitness s of two yeast strains. 




(19) 
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Figure 15. (Colour online) The radii of yeast colonies as a function of time. Yellow 
squares and black circles correspond to colonies of the wild-type and the advantageous 
sterile mutant, respectively. After an initial transient (time < 90h), described in more 
detail in the text, the radii arc well fitted by a straight line (red), in accordance with 
a constant expansion velocity. Inset: Instantaneous velocity ratio as function of time. 
The instantaneous velocities at a specific time arc determined from linear fits to the 
radii of the five surrounding time points. The black vertical line indicates the starting 
time for the fit in the main figure. 



5.1. Testing the equal-time argument 

We first tested the validity of tlie equal-time argument results of section H] with a 
particular pair of S. cerevisiae strains that have a large fitness difference: the wild-type 
and a faster growing mutant, which owes its advantage to the removal of a metabolically 
costly mating system; see the supplementary information (section SI) for more details. 
The equal-time argument assumes spatial expansion at a constant velocity. To see 
whether this assumption was valid in our experimental system, we measured the increase 
of the radius of circular yeast colonies over time, see figure [15l During the first day, 
colonies barely grow, presumably because it takes time for the populations to reach the 
carrying capacity in the spatial region of the inoculum. Afterwards, the colony fronts 
expand at about 20yum/h. The rate of expansion first slows down gradually over time, 
presumably because of nutrient depletion and drying out of the agar gel. For larger 
times, a stationary expansion front has established, and colonies grow at a constant 
expansion velocity. Since the latter is a pre-requisite for the equal-time argument to 
apply, we only consider these later times in our further analysis. A similar behaviour 
was observed for the expansion of yeast colonies from a linear inoculation (data not 
shown), although the absolute values of the velocities were very different; e.g. for the 
wild-type, we obtained V2 = 14.5 ± 1.9 /im/h for linear vs. V2 = 25.0 ± 0.5 fim/h for 
circular expansions. This difference is presumably due to larger colonies and therefore 
more severe nutrient depletion in linear inoculations. 

We next tested whether the equal-time argument correctly predicts shapes of 
boundaries between the strains: linear for linear expansions, logarithmic spirals for 
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Figure 16. (Colour online) Fitness estimation from linear expansion sectors, (a) 
S. cerevisiae colony grown from a linear inoculation at the bottom of the picture. A 
sector of the advantageous sterile mutant (black) emerges in the predominantly wild- 
type colony (yellow). The sector boundaries inferred from the image are shown with 
red lines. The scale bar is 500 //m. (b) Sector boundaries (blue dots) extracted from 
the image shown in (a), and fits (lines, > 0.995) to equation Note that the 
early part of sector growth differs from the later part. One possible explanation for this 
difference is the sector establishment process discussed in Sec.[3]when sector boundaries 
are not fully separated and interact with each other. The equal-time argument does 
not apply in this case. 



Figure 17. (Colour online) Fitness estimation from radial expansion sectors. The 
same as figure 1161 but for a circular geometry. In (a) , only the top half of a circular 
colony is shown. The smaller red circle shows the inoculum, and the larger red circle 
marks the colony radius. The scale bar is 1 mm. 



circular expansions, and circular for colliding colonies. This was indeed the case, see 
figures [16], [ITJ and [HI Note that in the case of sectors, the equal-time argument applies 
only at later times. During earlier times, the expansion velocity is not constant, as 
discussed above. In addition, during this initial transient, sectors have not yet fully 
established, see the discussion in section [3] and in the supplement section S5. Since both 
constant expansion velocity and fully established domains of each strains are required 
for the equal-time argument, we fit our theoretical predictions equations ((HI) and ([TOll 
to the sector shapes only for later times, finding excellent agreement (r^ > 0.995). 




0.5 
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Figure 18. (Colour online) Fitness estimation from colony collisions. The wild- 
type (yellow) colony meets the colony of the advantageous sterile mutant (black). The 
red lines are the fits of colony boundaries by circles. The relative fitness of the colonies 
can be measured from the radius and center of the circle fitted to the interface between 
the colonies; see equations (|18p and ([T5)) . The scale bar is 1 mm. 

5.2. Measuring relative fitness 

This agreement allows us to use the analytical results from the equal-time argument in 
order to measure the relative fitness s. From our analysis, the relative fitness can be 
estimated by five different methods, using (i-ii) ratio of expansion velocities vi and V2 
of isolated colonies in linear and circular geometries, (iii-iv) sector shapes in linear and 
circular geometries using equations (IE]) and ( fTOl) . and (v) the interface shape of colliding 
circular colonies using equations ( IT81IT91) . For consistency, we used the same time window 
for all assays, see the supplementary information (section SI) for details. To provide 
a reference, we also measured the relative growth rate during exponential phase in a 
well-mixed test tube, either in separate or in mixed cultures. 

The relative fitnesses s obtained from the different fitness assays yield quite similar 
results, summarized in table [TJ Indeed, the measurements are not significantly different 
from each other, except for the circular sector result which differs significantly from 
both the liquid culture competition and the colony collision assays {p < 0.05, see the 
supplement section S6 for details on the statistical testing procedure). The deviation for 
circular sectors could be caused by a systematic error in sector analysis. However, a likely 
explanation of the disagreement between different fitness estimates is some additional 
spatial structure not accounted for in our theory. Indeed, yeast colonies do not only 
expand on the surface of a Petri dish, but they also thicken over time to a height of 
about 1 mm, which is neglected in our two-dimensional theory. It is therefore possible 
that the advantageous mutant grows on top of the wild-type, producing an apparently 
larger sector and leading to an overestimate of Vi/v2. 

Table [U shows that the different fitness assays have standard deviations that vary 
over an order of magnitude, and therefore have very different accuracies. Expansion 
velocities of isolated single-strain colonies are the most straightforward measurement of 
fitness on a plate. However, they are less accurate than the sector and colony collision 
assays, as reflected by their high standard deviations in table [1] and the fluctuations of 
the instantaneous velocity ratio in the inset of figure [HI In sector and collision assays, 
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Method 


SplpptiArp ^1 H Ar^i n t ^1 cf p ^ 


Linear expansion 


velocity ratios 
sectors 


0.10 ±0.16 {N = 11) 
0.20 ±0.13 (iV = 23) 


Radial expansion 


velocity ratios 
sectors 

colony collisions 


0.16 ±0.08 (A^ = 19) 
0.23 ± 0.04 {N = 24) 
0.17 ±0.02 (A^ = 9) 


Liquid culture 


growth rate ratios 
competitions 


0.17 ±0.03 (iV = 3) 
0.18 ±0.02 (iV = 3) 



Table 1. Comparison of relative fitnesses measured by different methods. Errors are 
standard deviations (not standard errors of the mean) . The number of repHcates (iV) 
is given in parentheses. Note that the accuracy of different assays varies by about 
an order of magnitude. The large standard deviations for linear expansions arc, at 
least partially, due to front undulations, which make sector boundaries irregular and 
sector angles more variable. Sector and collision measurement have smaller standard 
deviations compared to direct velocity measurements. We attribute this distinction 
to the fact that both strains experience exactly the same local environment in the 
sector assay, but only approximately the same environment in the velocity assay. 
In particular, some environmental parameters-like the local dryness of the agar gel- 
are hard to control, and even identically prepared Petri dishes inevitably have slight 
differences in these parameters. Such variations affect velocity measurements, where 
the two strains are grown on two different (but identically prepared) Petri dishes, but 
do not affect sector measurements, where the strains are grown on the same Petri dish, 
and compete at the same point in space. 

the two competing strains are in exactly the same environment, and inevitable slight 
differences in the experimental conditions, such as humidity of the agar, influence both 
strains equally. This is not true for the isolated colonies of the expansion velocity assays, 
which are therefore more variable. Even more important, we found that relative fitnesses 
obtained from experiments on different batches of plates were significantly different when 
determined from expansion velocities, but not when determined from sectors or colony 
collisions; see the supplementary information (section SI). Nevertheless, the ratio of the 
expansion velocities is similar for linear and circular expansions, see table [Tj despite a 
large difference in the absolute values of the velocities, as discussed above. 

Linear expansion assays also have large standard deviations. This is probably due 
to undulations of linear fronts [26] , clearly visible in figure [H which distort the sector 
shapes. These undulations are significantly reduced in the circular geometry, leading to 
smaller standard deviations. It is thus advantageous to determine fitness from circular 
rather than linear expansions. 

From this discussion and table [H it follows that the radial expansion sector and 
colony collision assays are the most reliable assays to measure relative fitness in a spatial 
competition. We therefore only use these two spatial fitness assays in the following. 

We have so far found agreement between experiments and our equal-time argument 
predictions, but only for one pair of strains with one particular fitness difference. 



Selective sweeps in growing microbial colonies 



26 




Figure 19. (Colour online) Comparison of the selective advantage in well-mixed 
liquid culture and in spatial expansions. We varied the relative fitness using the 
drug cycloheximide for competitions of the cycloheximide-sensitive wild-type with a 
cycloheximide- resistant mutant. The fitness s measured with radial expansion sectors 
(red circles) and colony collisions (green squares) agrees well with the fitness Swm from 
the liquid competition assay, since all points lie close to the diagonal (black line, not a 
fit). The agreement of the spatial fitness s with the liquid fitness Swm is predicted by 
our theoretical model when migration is driven by cell growth, see the corresponding 
figure ini 

Therefore, we performed further experiments over a range of fitness values. To this 
purpose, we competed a strain resistant to with a strain sensitive to cycloheximide (a 
drug inhibiting translation) for varying cycloheximide concentrations in the medium. 
The relative fitness of the resistant strain, as e.g., measured with a liquid culture 
competition assay, increases linearly with the drug concentration, see supplementary 
figure S8. For all concentrations tested, radial sectors and colony collision boundaries 
could be well fitted with logarithmic spirals and circles, respectively, as predicted by 
our theory. More importantly, we were able to compare the fitness advantage in liquid 
culture to the fitness advantage on Petri dishes s over a wide range of relative 
fitnesses. We found good agreement between the two spatial assays and the liquid 
competition assay, see figure [191 

5.3. Testing the reaction- diffusion model 

Our reaction-diffusion model predicts an agreement of the spatial and well-mixed relative 
fitness, if migration is driven by cell growth, i.e. Dqi ~ gi and Dq2 ^ g2, as discussed at 
the end of section |3l The results shown in figure [19] are therefore a direct confirmation 
of this version of the reaction-diffusion model. In addition, the reaction-diffusion model 
predicts a constant front expansion velocity, given by equation [71 which is indeed 
observed in our experiments for large times, as shown in figure [151 Furthermore, the 
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reaction-diffusion model gives rise to the same, experimentally confirmed, macroscopic 
spatial pattern predicted by the equal-time argument, independent of microscopic details 
on the cellular length scale. We therefore preformed competition experiments with 

5. cerevisiae strains that have different cell division patterns, as well as with the 
bacterium Pseudomonas aeruginosa, see supplementary information (section S3). All 
experiments could be well described by our theory. 

In summary, there are experimental subtleties that our phenomenological theory 
does not take into account, such as the expansion velocity slowdown, yeast colony 
thickness, or possible strain interactions 0. All these effects could contribute to the 
slight differences of the fitness values determined by different methods, see table [H and 
figure [T9l Nevertheless, our theory describes the shapes of established sectors and colony 
collisions very accurately. Considering that the methods to determine relative fitness 
are very different, it is remarkable that the obtained fitness values are so similar, in 
particular results from well-mixed liquid culture and spatial growth on agar surfaces. 

6. Discussion 

Natural selection in well-mixed populations leads to selective sweeps of beneficial 
genotypes occurring exponentially fast in time. However, in spatially expanding 
populations, competition results in more complicated temporal and spatial patterns. 
Since both the advantageous and deleterious genotypes can spread into uncolonized 
territories, their competition can result in sectoring patterns like that shown in figure [61 
Sectoring spatial patterns provide an alternative fitness assay to the commonly used 
assays based on competition in a well-mixed environment of a test tube. This alternative 
facilitates spatial evolutionary experiments, which might contribute to understanding 
adaptations in different environments. 

The spatial assay may also be more accurate, provided front undulations, nutrient 
depletion, variations in agar wetness, and other experimental complications can be 
overcome. For small fitness differences, the assay could in principle acquire sensitivity 
because it measures ^Jv'l/vl — 1 = \/ s(2 -|- s) ~ yf2s instead oi V\jv2 = 1 + s [2^ . 
For large fitness differences, higher accuracy could also result from longer observation 
times, as the deleterious strain survives longer in spatial settings. A deleterious strain 
is eliminated linearly (linear geometry) or logarithmically (circular geometry) in time, 
unlike in the well-mixed environment, where it is eliminated exponentially fast. More 
important, a spatial competition experiment could be superior to a well-mixed one when 
used for screening for beneficial mutations. On a Petri dish, many beneficial mutations 

% The mutualistic or antagonistic interaction represented by the terms with ei and £2 in equation (j4|) 
could also change the relative fitness in experiments where the strains are in physical contact compared 
to experiments where the strains are grown in isolation. Three dimensional yeast colonies have a 
relatively large contact angle with the agar surface at the colony edge [35]. Therefore, the effect of such 
hypothetical interactions between the two strains might not be negligible if the density of yeast cells 
at the colony edge is not sufficiently small, as it would be if cell density decayed exponentially at the 
frontier. 
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can be assayed in parallel from expansions started by a small fraction of fluorescently- 
labeled mutagenized cells mixed with wild-type cells. In addition, each mutation is 
spatially isolated, and a dense aggregate of cells only a few generations away from the 
original mutation could be easily collected for future use. 

Our analysis of competition during range expansions has applications for 
evolutionary dynamics in spatially extended habitats as well. In particular, the 
predictions for one of the most important quantities in evolutionary dynamics, the 
duration of a selective sweep, is substantially different between spatial and nonspatial 
models. Deleterious genotypes persist much longer in spatial populations because the 
beneficial mutations spreading by Fisher waves may have to travel large distances 
needed, e.g., to engulf the wild-type populations, as in figure [T3l 
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